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Reacting fluids in the expanding Universe: A new mechanism for entropy production 
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It is shown that two reacting cosmological fluids, each of them perfect on its own, which exchange 
energy and momentum without preserving particle numbers, give rise to an entropy producing 
'reactive' bulk stress of the system as a whole, as soon as the detailed balance between decay 
and inverse decay processes is perturbed. This demonstrates explicitly that particle generation 
is dynamically equivalent to an effective bulk pressure. We derive a semiquantitative formula for 
the corresponding new kinetic coefBcient and evaluate the latter for the out-of-equilibrium decay of 
heavy, nonrelativistic particles into radiation. It turns out that the associated reactive bulk viscosity 
may be more than one order of magnitude larger than the conventional bulk viscosity, calculated, 
1^^ ' e.g., in radiative hydrodynamics. 
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(in ! I- INTRODUCTION 

Understanding the origin of the entropy content of the Universe is one of the major problems in cosmology. Entropy 

production requires deviations from thermodynamical equilibrium. Such deviations typically occur if microscopic 

interaction time scales become comparable vi^ith cosmologically relevant time scales. Phenomenologically, entropy 

producing processes may be described with the help of dissipative fluid models. From the standard dissipative fluid 

^-^ \ effects, heat conductivity, shear viscosity and bulk viscosity, the latter is the most favoured one in a cosmological 

r»o ■ context since it is the only effect which is compatible with the symmetry requirements of the homogeneous and 

f^ isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) universes. The existence of a nonvanishing bulk viscosity 

C^ for a simple relativistic gas has been known since the work of Israel ||l[ . A radiative bulk viscosity coefficient was 

^^ derived by Weinberg who also gave an analysis on its role in cosmology, followed by subsequent investigations 

(^ , by Straumann |g||, Schweizcr [^ and Udey and Israel ||5[]. Furthermore, bulk pressures were shown to play a role in 

O i' quark- hadron transitions 0. It was explicitly demonstrated recently [Q that the occurence of a cosmological bulk 

X \ viscosity is a general feature of any system of interacting fluids with different cooling rates in the expanding universe. 

^ ■ The considerations of all the mentioned papers were based on the assumption that the particle numbers of each of 

"^ ' the components are separately conserved during the interaction process. But particle number conservation is only a 

C^ . very special case, especially at high energies. In the present paper we consider a two-fluid universe with mutually 

L^ ' interacting and reacting components that are allowed to exchange energy and momentum and to convert into each 

• i-H , other. The corresponding change in the number of fluid particles is modelled by simple rate equations. In general, 

/^ ' these reactions do not preserve the total number of fluid particles. Each of the interacting and reacting fluids is 

H \ assumed to be perfect on its own. As long as the components are in thermal equilibrium and the conversion rates 

■ - - ' obey detailed balance relations, the entropy production of the system as a whole vanishes. Equilibrium situations like 

these are assumed to be realized at very high temperatures when the equations of state of both of the components 

and of the system as a whole are those for relativistic matter. As the Universe cools down, one of the subsystems, say 

fluid 1, becomes nonrelativistic while the other fluid, say fluid 2, continues to behave like radiation. The equations 

of state of the fluids are now different. Moreover, the detailed balance relations are no longer satisfied in general. 

The reactions may proceed preferentially in one direction, the inverse processes being suppressed. Typical cases for 

this kind of situation are the out-of-equilibrium decay of heavy bosons into quarks and leptons during the process 

of baryogenesis, or the scalar field decay into relativistic particles in the reheating phase of 'standard' infiationary 

scenarios g]. Also, the annihilation of electrons and positrons into photons at the nucleosynthesis energy scale belongs 

to this category. 

The purpose of this paper is to demonstrate that the difference in the equations of state of the fluid components 
and the simultaneous deviation from the detailed balance in the rate equations for the interfluid reactions give rise 
to a new entropy producing cosmological effect that manifests itself as a 'reactive' bulk pressure. We shall derive 
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here a semiquantitative formula for this effect for small deviations from equilibrium and apply it to a typical out-of- 
equilibrium situation in the early universe. 

To point out the difference to conventional bulk pressures, calculated, e.g., in radiative hydrodynamics, we emphasize 
that the physical reason for the effect under consideration here is twofold. The first reason is that because of their 
different equations of state the fluid subsystems are subject to different cooling rates in the expanding Universe. The 
second reason are the deviations from the equilibrium balances between the decay and inverse decay reactions of the 
fluids. It is the combination of both these causes that produces the effect to be discussed in the present paper. For 
conventional bulk viscosities, differences in the cooling rates alone are sufficient |^,0. 

The described effect implies that in a reacting mixture the generation of particles is dynamically equivalent to an 
effective viscous pressure. For a decay rate comparable with the expansion rate we show this effect to be more than 
one order of magnitude larger than the corresponding effect due to a conventional bulk viscosity. 

The structure of the paper is as follows: Section 2 establishes the basic theory of two interacting and reacting fluids 
with generally different temperatures. Assuming that the interactions establish an approximate equilibrium between 
the fluids, we investigate the same system within an effective one-temperature model in section 3. Comparing both 
descriptions, specifically the respective expressions for the entropy production, we derive a general, semiquantitative 
formula for the reactive bulk pressure of the system as a whole in section 4. Section 5 presents an application of this 
formula to a mixture of a nonrelativistic and a relativistic fluid and compares the reactive bulk viscosity coefficient 
with the conventional one. A brief summary of the paper and a few comments on cosmological applications are given 
in section 6. 
Units have been chosen so that c = ks = h = 1. 

II. THE TWO-FLUID MODEL 

The energy-momentum tensor T**^' of the cosmic substratum is assumed to split into two perfect fluid parts, 

rpik rpik I rj-iik / 1 \ 

with (A = 1,2) 

Tf = p,mV+p>^'=. (2) 

p^ is the energy density and p^ is the equilibrium pressure of species A. For simplicity we assume that both components 
share the same 4- velocity u\ The quantity /i*'^ is the projection tensor h'^'' = g''^ -I- u'u'^. The particle flow vector iV^ 
of species A is defined as 

K =n.u\ (3) 

where n^ is the particle number density. We are interested in situations where neither the particle numbers nor the 
energy- momenta of the components are separately conserved, i.e., conversion of particles and exchange of energy and 
momentum between the components are admitted. The balance laws for the particle numbers are 

K-. = ^A + en, - nj:, , (4) 

where 9 = m*j is the fluid expansion and F, is the rate of change of the number of particles of species A. There is 
particle production for F, > and particle decay for F, < 0, respectively. For F, = 0, we have separate particle 
number conservation. 

Interactions between the fluid components amount to the mutual exchange of energy and momentum. Consequently, 
there will be no local energy-momentum conservation for the subsystems separately. Only the energy-momentum 
tensor of the system as a whole is conserved. 

Denoting the loss- and source-terms in the separate balances by i^ , we write 

T% = -t\ , (5) 

implying 

P^+0(P^+Pj-"aC - (6) 

and 



{p,+P^Ju'^+P,^kh'^'^~htt\. (7) 

All the considerations to follow will be independent of the specific structure of the i^ . In other words, there are no 
limitations on the strength or the structure of the interaction. 
Each component is governed by a separate Gibbs equation 

r,ds, = d^+p,d— , (8) 

where s^ is the entropy per particle of species A. Using eqs.(^ and (^) one finds for the time behaviour of the entropy 
per particle 

^aTaSa == Uat°^ - {Pa +Pa)'^a ■ (9) 

With nonvanishing source terms in the balances for n^ and p^, the change in the entropy per particle is different 
from zero in general. Below we shall deal with the special case that the terms on the right-hand side of eq.(||) cancel. 
The equations of state of the fluid components are assumed to have the general form 

PA=PA{nA,T^) (10) 

and 

Pa ^ Pa i^A^T^) , (11) 

i.e., particle number densities n^ and temperatures T^ are regarded as the basic thcrmodynamical variables. The 
temperatures of the fluids are different in general. 

Differentiating relation (pi]), using the balances (0) and (0) as well as the general relation 

dp A _ Pa+Pa T^ dp^ 

that follows from the requirement that the entropy is a state function, we find the following expression for the 
temperature behaviour [p|-pT[: 

Both the source terms F^ and Uaf^ in the balances (Q) and (g) backreact on the temperature. The numerator of 
the second term on the right-hand side of eq.([l3|) coincides with the right-hand side of eq.(|9|), i.e., the corresponding 
terms disappear in the special case s^ = 0. 

For F^ = Uaf^ = and with 8 = 3R/R, where R is the scale factor of the Robertson- Walker metric, the equations 
of state Pr = UrkTr, Pr = SurkTr for radiation (subscript 'r') reproduce the well known T^ ^ R~^ behaviour. With 
Pm — rimkTm, Pm — rimm + ^UrnkTm One rccovers Tm ^ R^^ for matter (subscript 'ttt-'). 

The entropy flow vector S'° is defined by 

S: = UaSau'^ , (14) 

and the contribution of component A to the entropy production density becomes 

S''A,a^ri.^sj:^+n^k^ 

^aTa J T^ 

where relation (^) has been used. 

According to eq. (m the condition of energy-momentum conservation for the system as a whole, 

(Tf -f Tf ).^. = , (16) 

implies 

t° = -tl . (17) 



There is no corresponding condition, however, for the particle number balance as a whole. Defining the integral 
particle number density n as 

n = n^+n^ , (18) 

we have 

h + Qn = nV , (19) 

where 

nT = n,r, +n^r^ . (20) 

r is the rate by which the total particle number n changes. We do not require F to be zero since total particle number 
conservation is only a very special case, especially at high energies. 
The entropy per particle is |12| 



Pa +Pa Ma 



(21) 



where fi^ is the chemical potential of species A. Introducing the last expression into eg. ([Tsl) yields 

Sl.,a^-^n,T, + ^. (22) 

^ A ^ A 

For the total entropy production density 

SI = S^^, + 5^% (23) 

we obtain 

The condition 5".° = requires the well-known equilibrium conditions (see, e.g., |13| (chapter 5)) 

Ml = M2 > ^1 = ^2 . (25) 

as well as F = 0. 

From now on we assume that the source terms on the right-hand side of eq.(P) cancel among themselves, i.e., that 
the entropy per particle of each of the components is preserved. The particles decay or come into being with a fixed 
entropy s^. This 'isentropy' condition amounts to the assumption that the particles at any stage are amenable to a 
perfect fluid description. With s^ = 0, via eq.(P) equivalent to 

Uatl = ip^,+pjT,, (26) 

the expression dlS) for the temperature behaviour and the expression (p4h for the entropy production density simplify 
considerably. The condition (g^ establishes a relation between the source terms in the balances (^ and (j^) which 
originally are inde pen dent quantities. 
Combining eqs.([l7|) and (|2q), one has 

Uai° = {p,+p,)r, = -uJl ^-{p,+p,)T,, (27) 

which provides us with a relation between the rates Fj^ and F^ : 

T. Pi + Pit- 



P2 +P2 
Inserting the last relation into equation (Bfl) yields 



nV = rij^Fj^/ij 


■ 1 


1 " 



(28) 



(29) 



The quantities h^ = {p^ + Pa) /^a ^^^ the enthalpies per particle. Total particle number conservation, i.e., F = is 
only possible if the enthalpies per particle of both components coincide, i.e., for h-^^ = h^. 
With the relations (pq) and (E9h the entropy production density (04) becomes 



S?, = {p,+P,) 



Pi +Pi P2+ P2 



n, r, h 



1 '"1 



•^1 _ £2^ 



(30) 



We emphasize that acording to the equations of state (Jlfl) and (11) the quantities pi, p-^ and s^ depend on T^, while 
P2 , P2 and Sj depend on T^ . In general, we have T^ ^ T^ . 

From the expression (pG) it is obvious that a decay of fluid 1 into fluid 2 particles is accompanied by a nonvanishing 
entropy production unless s^/s^ = h^/h,^. For /i^ = /i^ = the latter condition reduces to T^ ^ T^. This situation 
corresponds to the case that all the entropy of fluid 1 is transferred into that of fluid 2. On this basis one usually 
discusses, e.g., the cosmological electron-positron annihilation at an energy scale of about 1 MeV ( [p|Jl^-|l^] . ) 

III. THE EFFECTIVE ONE-TEMPERATURE MODEL 

There exists a more familiar, alternative description of a two-component fluid close to equilibrium, that is based on 
a single Gibbs-equation for the system as a whole: 

TAs = d^ +pd- - (/i, - /i,)d^ , (31) 

n n ' n 

where s is the entropy per particle. The temperature T is the equilibrium temperature of the whole system. The 
temperatures T^ and T^ of the previous section do not appear as variables in the present effective one-temperature 
description. An (approximate) equilibrium for the entire system is assumed to be established through the interactions 
between the subsystems on the right-hand side of the balances (0) and (H). Furthermore, we assume that analogous 
to the equations of state (^0|) and { \\/\\ ) the cosmic fluid as a whole is characterized by equations of state 

p^p{n,n^,T) (32) 

and 

p = p{n,n^,T) , (33) 

where p is the equilibrium pressure and p is the energy density of the system as a whole. From the special case of a 
mixture of radiation and nonrelativistic matter discussed below eq.(|l^) it is obvious that the overall quantities p and 
p depend on two independent number densities which we have chosen to be n and n^ . As long as the pressures are 
those for classical gases, i.e., p^ = n^T, the equilibrium pressure p of the system as a whole depends on n = n^ -\- n^ 
only and the separate dependence on n^ on the right-hand side of eq.(p2) may be omitted. 



If in the expressions (21) for the entropies per particle the temperatures T^ and T^ are identified among themselves 
and with T, and p (T) = p^ (T) -I- p^ (T) as well as p (T) = p^ (T) + p^ (T) are used, the description based on relation 
( pi| ) is consistent with the description relying on the relations ^ for ns (T) — n^s^ (T) + n^s^ (T). 



The equilibrium temperature T is defined by 

Pi{ni,Tj+P2{n2,T.^)^ p{n,n^,T) . (34) 

As it was shown in 0] for F = F^ = F^ =0, this generally implies 

PAni,Tj+p^{n.^,Tjj^p{n,n^,T) . (35) 

For perfect fluids 1 and 2 with separately conserved particle numbers the difference between both sides of the latter 
inequality is a viscous pressure. Its existence is a consequence of the different temperature evolution laws of the 
subsystems. As it is obvious from eq. ([l3|) the coohng rate T^/T^ is different from T2/'^2 S'^sn for F^ = F^ = if the 
subsystems are governed by different equations of state. The expansion of the Universe tends to increase the difference 
between T\ and T^, i.e., to drive the system as a whole away from equilibrium. This manifests itself as a bulk viscous 
pressure [Q. 

We shall show here that nonvanishing source terms F^ in the particle number balances (0) give rise to a new type 
of effective bulk pressure, i.e., to enlarged entropy production. 



In order to find an explicit description of this phenomenon we shall use the following line of arguing. We anticipate 
that deviations from detailed balance, i.e., F^ 7^ 0, leading to F ^ in general, will generate an effective 'reactive' 
bulk pressure tt^^^^^ . For the corresponding energy-momentum tensor of the system as a whole we write from the 
outset 



T^ = pu'u" + ip + ^„„,, ) h"" . (36) 



To separate the reactive bulk pressure from any other dissipative phenomenon, we have ignored here the possibil- 
ity of conventional bulk pressures, dealt with elsewhere m and of nonvanishing heat fluxes and shear stresses in 

- ' ' ' - ' - — - — ik - - - - - 

■k 



inhomogeneous and anisotropic cosmological models. The relations T*^ — then imply the energy balance 



P + Qip + P + TT,,,J = (37) 

for the effective one-temperature description. Afterwards we shall determine tt^^^^^ by the consistency requirement of 
the one-temperature description on the basis of eqs.(|3l|), (0), dl9), and (p7|), with the two-temperature description 
based on eqs.(||), (Q), (||), ([l7|), and (26). More specifically, 7r^^„^t will be obtained by comparing eq.(23), leading to 



the two-temperature expression ( pO[) for the entropy production density with the corresponding expression for 5".° to 
be calculated below within the one-temperature picture. 

From the Gibbs-equation ( pi|) one finds for the change in the entropy per particle 

ns^-^n ^-^l^F-I^f^^^^i^VF^-Fj . (38) 

Even for s^ = s^ — we have s 7^ in general. From S*" — nsu"' the expression for the entropy production density is 

S^^ = nsT + ns . (39) 

Introducing here eg. (|38|) , we find 

S^. = 4-.... - ^F - I^ (t^^^) [F, - FJ , (40) 



T '''"* T n \ T 

where 



^^P±P.Ts (41) 



is the effective one-component chemical potential. 

Applying in eq.(En) the relation (ESh as well as the decomposition 



nfi = n^p.^+n^fi^ (42) 

for the chemical potential, the entropy production density may be written as 

,Qi rji 111 I / rji react \ / 

T \Pi+Pi P2+P2J T 

Here we used p{T) = p^ (T) -I- p^ (T) and p{T) = p^ (T) + p^ (T), keeping in mind that all the quantities, including 
those of the subsystems, depend on T, not on T^ . There is only one temperature in this kind of description, namely 
T. 

The expression (^3|) for the entropy production density is a sum of terms each being a product of thermodynamical 
'fluxes' and thermodynamical 'forces'. According to the standard phenomenological theory, the requirement S°^ > 
may be fulfilled by linear relations between 'fluxes' and 'forces'. The new aspect here, compared with the standard 
one is, that one has two scalar fluxes, namely tt^^^^j and F^. This will generally result in cross effects between both 
phenomena. A related discussion of chemical reactions in the expanding Universe was given in |lq ] . 

Within the causal thermodynamics of irreversible processes, both fluxes become dynamical degrees of freedom on 
their own and one arrives at coupled evolution equations for 7r,.^^^j and F^ . Investigations along these lines have been 
carried out for a one-component system by Gariel and Le Denmat |17| . 

In the present paper we do not follow this kind of arguing of the phenomenological theory. Instead, we shall find 
an explicit expression for 7r,.^^^j in terms of F^ with the help of semiquantitative arguments. The quantity connecting 
■""react ^^"^ ^1 rcprescnts a new kinetic coefficient, corresponding to the mentioned cross effect. 



Strictly speaking, the computation of tt^^^^j has to come from kinetic theory, which, however, is beyond the scope of 
the present paper. Our derivation may be regarded as intermediate between the phenomenological and kinetic levels 
of description. 

From eq.(p3|) the alternative expression for the entropy production density S""^ within the two-temperature model 
becomes 



^■a — "-1*1^1 



^1^1'' 1 I 



(44) 



with Sj = Sj (T^) and s^ = s^ (T^) from (|2l|), where we assumed each of the subfluids to be perfect on their own, 
i.e., S-^ = s^ ~ 0. The two-temperature formula (44) has to be consistent with the one-temperature formula (H3). 
This consistency requirement will be used to determine tt^^^^^ . The obvious difficulty one encounters here is that one 
has to deal with three generally different temperatures: The temperatures T^ and T^ of the components 1 and 2 in 
formula (B3) and the equilibrium temperature T of the system as a whole in (Bj). The dynamics of T^ and T^ are 
given by eqs.(|l^) with eqs.(|6|). In order to compare the expressions (^) and ([44|) for the entropy production we have 
to find a corresponding law for T. Such a law may be obtained via similar steps that led us to the expression (|l3|). 
There exists, however, the following complication. Because of the additional dependence of p on rij^ one has now three 
partial derivatives of p: (dp/dT)^ ^ , {dp/dn)j, ^ , and (dp/dn-^)^ j,. The requirement that s is a state function now 

leads to ( |l,|l9|) 



dp p + P T dp n^ 
dn n n dT n 

generalizing relation ([l^), and the additional relation 



d 



dp rr. d 



(45) 



(46) 



Using the Gibbs-Duhem relations 



<^Pa ^ri^s^dT^+n^dn^ 



(47) 



for T^ — T together with eq.(|21|), one finds 



^^ dT n. 



(48) 



Consequently, the relations (Ba) and (|6h may be written as 



dp 


P + P 


T dp 


n. 


' Pi 


P2 


dn 


n 


ndT 


n 


.^1 


""■2 



(49) 



and 



dp_ 

dn. 



Pj_ _ fh 
n, ri„ 



(50) 



with Pj = Pj (T) and p, — p^ (T), since we are within the one-temperature description. 
Differentiating eq.(p3h, using eqs.(H), (n9h, (37), (Es), and (50), we obtain 



^T = - (6 - r) T^ - [Qn^, + r (p + p)] 

n, n„ J 



(51) 



Applying here eq.(|28|), following from s-^ — s^ = 0, and restricting ourselves from now on to classical fluids with 
p^ — n^T, we find the following evolution law for the equilibrium temperature: 



T 
T 



(e-r) 



dp Ott^, 



dp Tdp/dT ' 



(52) 



where the abbreviation 

dp _ dp/dT 
Wp ^ dp/dT 



(53) 



was used. 

The temperature law (p3) is a new result. It describes the behaviour of the equilibrium temperature of an interacting 
and reacting two-fluid mixture with nonconserved total particle number. It is different in its structure from the laws 
(O) with (Eq) because of the tt^^^^^ term. The resul t (p^ ) differs also from previous one-component approaches to 
particle production in the expanding Universe |10| , |l l| , |20|| 7^ The point here is that the 'isentropy' conditions s^ — 0, 
equivalent to relation (Eq) of the present paper are physically different from the case of 'adiabatic' particle production 
|10| , |ll| , P0| , characterized by s = 0. We conclude that 'adiabatic' particle production in the sense of 0,|l 20 



IS not 



consistent with the requirement of a perfect fluid behaviour of the subsystems in a reacting two-component mixture. 

IV. THE EFFECTIVE VISCOUS PRESSURE 



With the formulae (13^ and (52) we know the evolution laws for all three temperatures. We are now prepared to 
compare the expressions (Hs) and (H4) for the entropy production densities. To that purpose we write s^ — s^ (T^) 
and s^ = Sj (TJ as {A = 1, 2) 

s, (TJ = ,s, (T) + As, . (54) 

Using this splitting in eq. ( [44[) and comparing with eq. ( [43| ) yields 

- ^TT^^^ = r^n^ As^ -I- r,n,As, . (55) 

Since according to eqs.(O) and (B3) the temperatures T, and T are different, the quantities As, = s, (T, ) — s, (T) 
will be nonvanishing in general. However, as long as the detailed balance relations F, = are fulfilled, 7r^_,^^j vanishes 
even for As, ^ 0. (We recall that this case leads to conventional bulk stresses, discussed elsewhere 0). 

In order to find an explicit expression for tt^^^^^ we use the following heuristic mean free time argument (cf. [[7|). 
Let T be the characteristic mean free time for the interaction between both components, generally symbolized by i" 
in eq.(H). The time r is assumed to be much larger than any characteristic interaction time within each of the fluids. 
Consequently, the latter may be regarded as perfect on time scales of the order of r. The interaction between the 
fluids is pictured by 'coUisional' events, where r plays the role of a mean free time between these 'collisions'. During 
the time interval t, i.e., between subsequent interfluid interaction events, both components then evolve according to 
their internal perfect fluid dynamics, given by eqs.(0), (@) and (O) with condition (Eq). 

Let us now assume that at t = i,, there exist deviations from detailed balance due to the circumstance that 
one of the components, say fluid 1, becomes nonrelativistic and starts to decay into the relativistic fluid 2. Let us 
further assume that both components were in thermal equilibrium aX t = t^, due to their mutual interaction, i.e., 
Tj (i^) = Tj (t^) = T{t^), equivalent to As^ (tg) = As^ (i,J = tt^^^^ (i,,) = 0. Until a subsequent interaction event 
after a mean free time r between the 'collisions' the components evolve effectively free due to their perfect fluid 
dynamics. According to the temperature laws (O) and (p2|), the temperatures T^, T^ and T evolve differently for 
t > t^^. At the time t with t — t^ < t we will have T y^ T, in general with T^ ^ T^. Consequently, the quantities Asa 
will no longer vanish at the time t. Nonvanishing values of As,, however, together with F, ^ 0, will give rise to an 
effective bulk pressure according to eq.(|55|). 

Up to first order in the temperature deviations we may write 

As, = (T,-r)('^) +... . (56) 



Applying the general formula 



ds._ 1 dp, ^^^^ 



dT n,T, 9T, 



that follows from the Gibbs-Duhem relation (p7|), we get 



- ^K.... = n,r,cl'' {T, -T)+ nj:J^' (T, - T) , (58) 



where 



n.. dT 



(59) 



are the specific heats of the components. Introducing the relation (ESh between the rates F^ and F^, eq.(p8|) may be 
written as 



©""^cct =r^ {p^ +pj 



n,Cv {T^-T) n^c^ (T^-T) 



Pi +Pi 



P2 +P2 



rijFj/ij 



c:'iT,~T) ^J{T,-T) 



(60) 



We recall that this reactive bulk pressure is a cross effect between the deviation from detailed balance, represented 
by Fj 7^ 0, and the deviation from thermal equilibrium, i.e., T^^T^^ T. The conventional bulk pressure is only due 
to the different cooling rates. Here, both deviations from the detailed balance and from the thermal equilibrium are 
necessary. 

Let again component 1 be the decaying component (i.e., F^ < 0) with an equation of state close to that for matter. 
Because of dp^/dp^ > dp/ dp > dp^/d£^, the temperature T^ cools off faster than the equilibrium temperature T, 
while T^ decreases slower than T [cf . ([l3|) , (|5^)]. Consequently, we expect T^ < T and T^ > T. In total, the result 
according to eq.(60) will be a negative bulk pressure corresponding to a positive effective bulk viscosity, in accordance 
with the second law of thermodynamics. 

Restricting ourselves to small temperature differences we have 



T,=TAt„) + it-t„)t,+ 



(61) 



and 



r = T(ij + (t-ijr + 



(62) 



Here, 



T^ and T have to be evaluated at t 



t,-,. While T^ may be taken from eq.(L3[) immediately, we have still 
to show that T (t^) is consistent with tt^^^^j (tg) = 0. This may be done by applying the relation (p3) that defines 
the equilibrium temperature T a,t t > t^, keeping in mind that the interaction between the fluids, generally given by 
the source terms in the balances (^ and (|7[), is modelled by 'collisional' events, with a mean free time r between 
the 'collisions'. During the interval (t — tg) < r, i.e., between subsequent interaction events, both cornponents are 
effectively free, i.e., governed by their internal perfect fluid dynamics according to eqs.(P, (g), (|l^), and (p6|). Assume 
again that through the interaction an element of the cosmic fluid is in equilibrium at a proper time tg at a temperature 
T{tg) = T^ (tg) = T^ (tg) withp{tg) = p^ (tg) +^3 ^t gVnev B , p{tg) aud p^ {tg) are shorts for p[n{tg) ,T (tg)] and 
Pa ["-A (^0) ' ^-4 (^0)]' respectively. Using the condition (p3) at the proper time t up to first order in t — tg, i.e., with 



PAt)^PAitg) + {t^tg)p,{tg) + .. 



(63) 



and 



Pit)^p{tg) + it-tg)pit0) + - 



(64) 



where p^ (tg) = p^ [n^ {tg),T^(tg}\ and p(ij = p[n{tg),n^ (ij,r(ij], applying eq.(|) on the left-hand side of 
eq.(|33) and using the relations (Es) and (pO|) on its righ-hand side, one finds 



T(ij = -r(e-F) 



dp 
dp 



(65) 



where all the quantities on the right-hand side of the last equation have to be taken at the point t — tg. Comparing 
with the temperature law (p2|), this proves the consistency of our initial configuration with it t (^o) — ^'^^ ^i (^o) = 
T, itg) = T {tg), i.e.. As, = [cf. eq.(|p)] and F, (ij ^ 0. 

Applying now the evolution laws (|l3|) and (|6^) at t = tg we get for the first-order temperature differences 



T 



-eT{t-tg) 



T, 



dpj 
QJ dp. 



T\ dp 



ej dp 



(66) 



The first-order approximation is valid as long as the conditions 
\ld — r ^\ {t — tg) < 1 are fulfilled. Recalling that an effective one- fluid description requires r < H^^, where H — 8/3 
is the Hubble parameter, one may discuss three different cases. 

(i) |r^ |~^ < r < H~^. There is a considerable change of the particle numbers on time scales of the order of the mean 
free time between the intcrfluid interaction events. 

(ii) T < H^^ < |r^|^^. Here, the characteristic times |r^|~^ for the particle number changes are larger than the 
Hubble time. 

(iii) T < Ir^l^^ < H^^. This case is intermediate between (i) and (ii). The particle numbers change substantially 
during one Hubble time but are approximately constant on scales of the order of r. 

Inserting the temperature differences ( pq ) into eq. ( |60| ) , the reactive viscous pressure may be written as 



, r rr.dp^ 



Qjdp, \ ejdp 
r^\ ap^ _/ r \ dp 
Qjdp, \ ej dp 



(67) 



where rj = t — t^. It can now be seen explicitly that 7r^_,^^j is always negative. For 'ordinary' matter dp^/dp^ lies in 
the range 1/3 < dp^/dp^ < 2/3. The lower limit corresponds to radiation, the upper one to matter, dp/ dp will take 
a value intermediate between dp-^/dp^ and dp^/dp^. The typical case is that massive particles decay into relativistic 
ones. Identifying again fluid 1 with the massive decaying component, we have F^ < 0. Furthermore, dp/ dp < dp^/dp^ 
and dp/ dp > dp^/dp^. Because oi dp^/dp^ > dp/ dp one has (1 — Fj/8) dp^/dp^ > (1 — r/6) dp/ dp, since a negative 
Fj enlarges the left-hand side of this i neq uality and a positive F diminishes its right-hand side In the second term one 
has dp^/dp^ < dp/ dp. The relations ( [28|) and ( P9| ) may be combined into 



ri2_ 
n 



1 



K 



(68) 



One has F < F^ for any F^ > 0. It follows that the factor (1 — F^/O) reduces the dp^/dp^ term more than (1 — F/0) 
reduces the dp/dp term. Consequently, (1 — T^/Q)dp^/dp^ < (1 — T/Q)dp/dp. Both contributions to tt^^^^j in the 
expression (|6^) are negative, i.e., the reactive bulk viscosity coefficient is positive, as is required by the second law 
of thermodynamics. Inserting the expressions (Eq) and (59) into equation (CT), the reactive viscous pressure may be 
written as 



T -^rjU-^h-^T 



-^.r,-i<hlT 



+ 



dp, 
dp, 

(1) 



dp 
dp 



dp2 
dp2 



dp 
dp 



dp, _ r\^ f I l_\ dp 

h^ h.^ J dp 



iK \dp, n 



n^hl 



dP2 



1 



1 



dp 

d~p 



(69) 



This is the general formula for the reactive bulk pressure in a two-component cosmological fiuid for small deviations 
from thermal equilibrium. It is the main result of the paper. By eq.(p39|) the quantity tt^^^^j is given in terms of F^ 
or, via eq.(|29|), in terms of the production rate F. This demonstrates explicitly that particle production in a two- 
component mixture is equivalent to an effective bulk pressure. The existence of a formula like eq.(p3) may be regarded 
as a semiquantitative justification for using effective viscous pressures in modelling particle creation processes on a 
phenomenological level. We recall that different from previous studies ( P, [lO|j2C| -|2^ ) the production process here is 
nonadiabatic. 

In the following we shall apply formula (p9h to the out-of-equilibrium decay of a nonrelativistic fluid into a relativistic 
one. 



V. THE OUT-OF-EQUILIBRIUM DECAY OF HEAVY PARTICLES 

Decay processes of heavy particles are supposed to play a role, e.g., during baryogenesis where hypothetical heavy 
bosons decay into quarks and leptons. Another application is the scalar field decay into relativistic particles in the 
reheating phase of standard inflationary scenarios (see, e.g., M). 
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We assume fluid 1 to be described by the equations of state for nonrelativistic matter, i.e., 

3 



while fluid 2 is a classical relativistic gas: 



The corresponding enthalpies are 



P2 = Sn^T^ , p^ = n^T^ 



(70) 



(71) 



and the specific heats 



K=m+ -T, 



(1) 3 



K = 4r, , 



c„ = 3 



(72) 



(73) 



The expressions for dp^/dp^ and dp/ dp become 



dp. 


2 


dp2 


1 


dp 


2 n, +71^ 


dp. 


3 ' 


dp2 


3 ' 


dp 


3n, + 2n^ 



Inserting eqs.d70)-d7l) into eq.(p3) and taking into account m ^ T we obtain 

nim n. 



4 n, + 2n„ 



1 



n,m r. 



An^T e 



T,r] 



Since p + p — p, + p^ + p, + P2 ^ n,m + -in^T, the ratio tt^^^^J {p + p) is 



1 n,m 



p + p An^rn + ^n^T n, + 2n2 



n,m Ir, 



An^T e 



Ti?? 



(74) 



(75) 



(76) 



In order to appropriately apply this formula we have to take into account the limitations following from the linear 
order restrictions (pi]) and ( |62[ ) for the temperature evolution laws. From eq.dl3) with ( p7| ) we find the linear relations 
( pi] ) and (|6^) to make sense for 



|e-rj(i-ij<i 



(77) 



To be cosmologically relevant, we expect JFJ to be of the order of 8, i.e., —T, = JFJ = a<d, where a is a positive 
constant of the order 1. For component 1 the linear evolution law is valid up to time intervals 



(t-h) 



1 



e+|Fj a + V '' 



(78) 



Assuming furthermore the nonrelativistic component 1 to dominate initially, it follows from eq.(28) that F^ > |F^|, 
such that Fj > is possible. In the latter case combination of ^Tq) with (U7t) for A — 2 yields 



a + 2, 



(79) 



Comparing the last expression with cq.(E8^ we get, by using the equations of state (fTG) and (|7l|). 



Pi 



4a + 2 
3 a 



-P2 



(80) 



The limitations of the linear approximation allow the energy density of the matter component to be larger than that 
of the radiation component by the factor 4 (a + 2) / {3a). For awl the ratio P1/P2 ~ 4 is allowed, while for a w 1/4 
we have p^/p^ ~ 12. Since eq.(pO[) is equivalent to 



n,m a + 2 
An^T a 



(81) 
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we find from eq.([75|) 

n, TO n, 



4 n, + 2n. 



(a + 3) a7]Q , (82) 



corresponding to an effective bulk viscosity coefficient ^.^^^^ 

C_=.^;^^^a(a + 3). (83) 

It is interesting to compare this bulk viscosity coefficient with the conventional one, characterizing the viscous pressure 
in a two-fluid system with conserved particle numbers. In general, the latter is given by R] 

^— 1^(1: -I) (I:- 1) ■ <-' 

which, for matter dominance, coincides in good approximation with Weinberg's expression for the bulk viscosity in 
radiative hydrodynamics [||j^ . With the equations of state (|7y) and (^l|) the expression ( |84| ) reduces to 

3 rij + 2^2 
Assuming now r]K, t (see the discussion below eq.(pO)) we find 

^«3£g.a(a + 3), (86) 

and, with eq.(|8l|). 



c 



3 (a + 2) (a + 3) . (87) 



For a « 1 one has Ceoct ~ 36C, while for a w 1/4 the quantity Crcact i^ more than 20 times larger than C^. The 
reactive bulk viscosity exceeds the conventional one by more than one order of magnitude and, consequently, provides 
the dominating contribution to the bulk viscous part of the entropy production during the decay process. 

VI. CONCLUSIONS AND OUTLOOK 

Particle number nonconserving reactions between two cosmological fluids, each of them perfect on its own, neces- 
sarily imply the existence of entropy producing reactive bulk stresses. This result of the present paper provides a 
semiquantitative justification for the frequently used approach of regarding particle creation processes as phenomeno- 
logically equivalent to effective viscous pressures. Different from most of the previous applications, the particle 
production turned out to be 'nonadiabatic' in the present case, i.e., the entropy per particle of the cosmic fluid as 
a whole does not remain constant. For small deviations from thermal equilibrium we derived a general formula for 
the corresponding new kinetic coefficient. The latter was explicitly evaluated for the out-of-equilibrium decay of 
nonrelativistic fiuid particles into relativistic ones. The reactive bulk viscosity coefficient in this case was shown to 
be larger than the conventional bulk viscosity by more than one order of magnitude. 

Our results indicate that bulk pressures are a general cosmological phenomenon. While not taken into account in the 
standard perfect fluid approaches, one may expect them to have considerably influenced the evolution of the early 
universe. Of special interest in this respect are far-from-equilibrium situations. The existence of a new type of bulk 
pressure in addition to the 'conventional' one may shed new light, e.g., on the question whether bulk pressures via 
their backreaction on the cosmological dynamics might give rise to a phase of inflationary growth of the cosmic scale 
factor (see p4| , [ll| , p5| and references therein). 

An out-of-equilibrium period that recently has attracted much attention (see, e.g. , |2^,|2^) is the reheating phase 
of inflationary universe models. On the basis of the temperature evolution law (|5^) it was shown to be possible to 
interrelate the description of nonequilibrium processes within the causal, second order Miiller-Israel-Stewart theory 
with particle creation during the 'preheating' stage pq]. A further subject of interest in this connection may be the 
'reheating' from primordial black hole evaporation [p9|. 
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